GRS 1915+105 as a Galactic Analog of a Fanaroff-Riley II Quasar 
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ABSTRACT 

We study the long term time averaged kinetic luminosity, Q, of the major 
flares of the Galactic microquasar GRS 1915+105 and the relationship to the 
intrinsic X-ray (bolometric) luminosity, Lbol, and scale it to that of a complete 
sample of SDSS/FIRST FR II quasars. If the scale invariance hypothesis for 



Jl _ 

black holes (BHs) holds then we show that the expected distribution in the Q - 

L bol scatter plane of GRS 1915+105 is consistent with FR II quasars for distances 

D = 10.7 - 11 kpc. We compare the specific values of kinetic luminosity and L bol 

during flares of GRS 1915+105 to that predicted by several 3-D MHD simulations 

of BH accretion flows with relativistic ejections. If FR II quasars are a scaled 

up version of GRS 1915+105, the data are consistent with numerical models 

when they contain an ergospheric disk jet and the BH spin is a/M = 0.99 or 

^ | a/M = 0.998 (we estimate a/M > 0.984). In the framework of scale invariance 

of BHs, our results may imply that FR II quasars also hold rapidly rotating BHs. 

O 

Subject headings: Black hole physics — magnetohydrodynamics (MHD) — galax- 
ies: jets — galaxies: active — accretion, accretion disks 

1. Introduction 

The black hole candidate GRS 1915+105 is well known for launching superluminal 
radio flares out to large distances at a much larger rate than any other Galactic object 



( JMirabel and Rodriguezlll994t iFender et alHl999t iDhawan et alll2000l ). The X-ray luminosity 
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of GRS 1915+10 5 is one of the highest of any known Galactic black hole (BH) candidate 



( jDone et alll2004l ). The existence of both relativistic outflows and high continuum luminosity 
make it tempting to speculate that GRS 1915+105 might be a scaled down version of a radio 
loud quasar. This is of particular importance because the time scales for radio evolution 
are reduced from the AGN (active galactic nuclei) time scales by many orders of magnitude. 
Thus, unlike quasars, it is in principle possible to see the details of the connection between the 
putative accretion flow (X-ray luminosity) and the superluminal jet launching mechanism. 

In this article, the idea of GRS 1915+105 as a scaled down FR(Fanaroff-Riley) II quasar 
is explored. The long term time averaged power of the relativistic major flares in GRS 
1915+105, Q, and the luminosity associated with viscous dissipation in the accretion flow, 
Lboi, are re-scaled in order to compare and contrast with the distribution of Q and Lboi of a 
complete sample of FR II quasars. These efforts are rendered credible by the recent study 
of the power required to launch individual major flares and their accretion state (intrinsic 
X-ray luminosity or bolometric luminosity, Lboi) just hours and minutes before ejection and 
during the brief 1 to 7 hour ejection event (Punsly and Rodrgiuez (2013), PR13 hereafter). 
The relevant results from PR13 that are required to perform this re-scaling are indicated 
in Section 2. In Section 3, this is compared to FR II quasars. In Section 4, the results 
of Sections 2 and 3 are considered in the context of numerical simulations of 3-D MHD 
(magnetohydrodynamic) accretion around spinning BHs. 



2. Estimating the Energy Output from Relativistic Ejecta in GRS 1915+105 



It was determined in iPunslyl (120121 ). P12, that knowledge of the time evolution of the 
spectral shape associated with a changing synchrotron - self absorbed (SSA) opacity, r, 
greatly enhances the accuracy of plasmoid energy estimates (constrains the size). The fre- 
quency and the width of the spectral peak provide two added pieces of information at each 
epoch of observation beyond the single epoch spectral index and flux density of the opti- 
cally thin high frequency tail that is traditionally used to estimate the ejected plasmoid 
energy. The evolving r combined with baryon number conservation, energy conservation, 
synchrotron cooling times and X-ray luminosity were used in P12 to eliminate uncertainty 
in the energy estimates. Namely, the proton content is minimal and a near minimum energy 
condition, E min pa m e c 2 , is shown to occur when the optically thin flux at 2.3 GHz, Sthin(2.3), 
is near maximum and r pa 0.1. In PR13 we assume that the detailed modeling of the time 
evolution of the flares from P12 can be used as a template for the time evolution of other 
plasmoids with less supporting data. This determines a set of equations for each flare in 
PR13 that are solved numerically with 4 inputs, the peak <Sthin(2.3), the spectral index of 
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the optically thin emission, D, and the Doppler factor, S. Figure 1 is a plot of the plasmoid 
energy estimated in Table 2 of PR13 as a function of the estimated peak £thin(2.3) for the 
major flares that were listed in Table 1 of PR13. Except when noted, a fiducial distance 
to GRS 1915+105 of D =11 kpc is assumed throughout the manuscript. Unlike Table 2 
of PR13, the energy is measured in the observer's reference frame instead of the plasmoid 
reference frame. The power law fit in Figure 1 w as made with the method of weighted least 



squares with errors in both variables (JReedlll989l ). The corresponding power law fit for each 



D value is used here to estimate flare energy from a single input, the peak 5thi n (2.3). The 
approximately daily 2.3 GHz and 8.3 GHz monitoring with GBI (Green Bank Interferom- 
eter) provides a database from 1996 through 2000 of 1967 days (taking account of gaps in 
coverage) for which one can determine Sthm(2.3). The resulting distributions of major flare 
energy and number from 1996 through 2000 are plotted in Figure 2. 

There are three major sources of uncertainty in the long term cumulative energy, E, 
from the flares 

1. The uncertainty in the estimate of the plasmoid energy 

2. The minimal Sthm(2.3) that is indicative of a superluminal ejection 

3. Possible short (weak) flares missed in the intra-day gaps in coverage 

The stochastic error in the first item is ignored because there are sufficient flares to 
drive the propagated random error to approximately 0. We consider every plausible value of 
D and the corresponding dependent 5 because of the large systematic uncertainty in 5. In 



order to estimate S from D, we assume that the kinematic results from iFender et all (119991 ) 



are common to the entire time frame from 1997 to 2000 as evidenced by interferometric 



observations of multiple flares (jDhawan et all |2000| ; iMiller- Jones et all 120051 ). The intrinsic 



spectral luminosity is S , t hi n (2.3)(5~^ 3+a ^. As D is varied from 10.5 kpc to near the maximum 
kinematically allowed value of 11 kpc, the intrinsic spectral luminosity changes by a factor 
~ (0.31/0.54) 3 - 9 = (1/8.7) which equates to a reduction of the plasmoid energy by a factor 
of 5 - 6. Alternatively, the discussions to follow can be phrased in terms of 5 instead of D. 
For item 2, we assume that the weakest flares have a peak Sth™ (2.3) = 30 mJy since this is 
the smallest value that can be clearly discerned from a background consisting of core flux 
variations and previous fading flares. Extrapolating the distributions in Figure 2 to lower 
cutoffs shows that the total energy output is rather insensitive to the low energy cutoff of the 
flares - most of the flares are weak, but they carry a small fraction of the total energy output. 
We estimate the contribution from weak flares in the intra-day gaps and the uncertainty from 
points 2 and 3 to compute Q, 
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Ejected Energy versus Peak Optically Thin Flux (D = 11 kpc) 



E = 4.95 x 10 38 [S thin (2.3)/1 mjy] 1 408 ergs 
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Fig. 1. — This scatter plot of the estimated energy of superluminal ejections from PR13 
evaluated in the observer reference frame versus the estimated peak Sthm(2.3) 
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Q 



cumulative energy of major ejections 
1967 days 
1.78 ± 0.21 x 10 45 ergs 



1967 days 



1.04 ±0.13 x lO^ergs/s 



3. Scaling to a Supermassive Black Hole 



In this section, GRS 1915+1 05 is re-scaled to FR II qu asar values of Lboi- A complete 
sample of quasars was created in iPunsly and Zhang I (1201 ll ) using the combined SDSS and 
FIRST databases. The radio sensitivity was adequate to detect extended emission below the 
FR I/FR II divide. Thus, a complete distribution of FR II quasars is attained. Furthermore, 
the Q estimates are far more accurate than other treatments of FIRST data in the literature. 
Radio images were used to subtract the jet emission and core emission on scales less than 
20 kpc in order to more accuratel y determine the op tically thin lobe flux density which is 
the most robust estimator for Q (jWillott et al.lll999l ). The blue dots in F i gure 3 are the 
Q-Lboi scatter plot from the top left frame of Figure 3 in IPunsly and Zhang I (]201ll ) with the 
conversion from the integrated optical/UV continuum luminosity to Lboi given by Equation 
(3) of that paper. The Lboi estimate is based on flux that is emitted ~ 10 6 — 10 7 years after 
the preponderance of plasma that is responsible for Q was eject ed from the central engine, 
based on lobe separation and estimates of lobe advance speeds (jWillott et al.lll999l ). Thus, 
the epochs are so displaced in time that there need not be a causal connection. In order 
to compare GRS 1915+105 to the AGN data, one needs to compare Q in GRS 1915+105 
to Lboi at epochs that are not likely to be causally related, i.e. randomly selected from the 
historical distribution of L D oi- 

This random data sampling is explored with Monte Carlo simulations. First, we create a 
distribution of quasar Lboi from the SDSS/FIRST data scatter, /quasar (Lboi)- Second, the Lboi 
distribution of GRS 1915+105 from PR13 Figure 16 and Q from Equation (1) are combined 

log (#7 



to create a random distribution of Q/L^i, /grsi9±5 log ( jf— ) . Then we create pairs of 
points representing the re-scaled GRS 1915+105 based on re-scaling Lboi to the quasar 
level. This is accomplished by randomly generating L bo i from /quasar (Lboi) then creating Q 



randomly for each L bo i from / G rsi915 



log 



Lhoi 



. The red dots in Figure 3 represent Monte 



Carlo simulations of the re-scaled GRS 1915+105 with Lboi distributed similarly to the 
SDSS/FIRST FR II sample. The explicit expressions for the quasar luminosity distribution 
(the blue dots in Figure 3), /quasar (Lboi), is consistent with a log-normal distribution, Z, with 
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Fig. 2. — The distributions of number and energy ejected in major flares emitted from GRS 
1915+105 from 1996 to the end of 2000 (see text for details). The black (red) curve is the 
cumulative distribution of flare number (energy) in the 1967 days, 
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a mean logarithm (in units of erg/s) of 45.96 and a standard deviation, 0.29, 

/ quasar [log(L bol )] = Z{p = 45.96,a = 0.29) . 



(2) 



The distribution of Lboi for GRS 1915+105 is well described by a log-normal distribution in 
Figure 16 of PR13 



/GRSi9i 5 [log(L bol )] = Z(ji = 38.73, a = 0.15); , D = llkpc 
From Equations (1) and (3), log(Q) = -1.71 + /i[/ G RSi9i5[log(£boi)]], thus 



/. 



GRS1915 



log 



Q 



L 



bol 



Z(fi = -1.71, a = 0.15) , D = llkpc 



(3) 



(4) 



Similar expressions for D < 11 kpc follow from Table 2 and the methods of PR13 (see Figure 
2 also). 

For D = 11 kpc the 2-D distribution of the re-scaled GRS 1915+105 is clustered near 
the peak of the quasar distribution. The smaller dispersion of the simulated data is expected 
because all points are generated by the same central engine as opposed to a variety of central 
BH masses and enveloping environments that are responsible for the quasar generated scatter. 
At D = 10.7 kpc, < 1/2 of the simulated data is consistent with FR II quasars and at D= 
10.6 kpc the distributions have become very distinct. Systemat ic uncertainty is fou nd by 
comparing the high biased quasar Q estimates i n Figur e 3 fro m IWillott et al.l (119991 ) with 
the low biased Q estimates from the methods of iPunslyl (120051 ). that are based on different 
assumptions. This yields a systematic uncertainty in log Q of 0.38 ± 0.08, too small to affect 
the implications of the Monte Carlo simulations. There is very little systematic error in 
Lboi since it is derived from the integrated continuum near the peak of the spectral energy 
distribution. The systematic uncertainty in Q for GRS 1915+105 was discussed in Section 
2. The main systematic error in Lh n \ is in the column density, Nh, to the source. From 
Belloni et all (119971 ) ; iMuno et all (119991 ) we expect a systematic error less than a factor of 2 
arising from the uncertainty in Nh, too small to affect our conclusions. 



4. Results in the Context of Simulations of Black Hole Accretion 

We explore the main results from PR13 in the context of numerical simulations 



1. Strong flares are launched when Lboi is at an elevated level (sometimes approaching 
the Eddington limit). 



FRII JET POWER AS A FUNCTION OF BOLOMETRIC LUMINOSITY 
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Fig. 3. — Monte Carlo simulations are used to compare the 2-D distribution in the Q - L\ iQ \ 
scatter plane for a complete sample of SDSS/FIRST FR II quasars (blue) and the re-scaled 
GRS 1915+105 (red) 
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2. During the 1 to 7 hours of major flare ejections, the time averaged power, < Q >, and 
time averaged intrinsic radiative luminosity, < Lboi >, are highly correlated. 



Our estimates of Lboi are based on models where the main contribution is due to thermal 
Comptonisation of soft (cold ~ 0.2 keV) photons by hot (~ 20-100 keV) electrons present 
in a so-called corona. The important parameters to estimate Lboi are kT in j, kT e , r and the 
Comptonised normalization (see Section 4.2.2 of PR13). < Q > / < Lboi > f° r ea ch D is 
computed from the values in Tables 2 and 3 of PR13 for each flare. The solid line in Figure 
4 is < Q > J < Lboi > averaged over all the flares for each D. The dashed lines represent the 
standard deviation, ±<x. The small dispersion argues strongly for an approximately constant 
ratio < Q > / < Lboi > at each D. We assume that this is the case and the errors in the 
individual < Q > and < Lboi > are artifacts of imperfect data and our estimation methods. 
Equivalently, the uncertainty in < Q > / < Lboi > for each D is a. 

There are two topologically distinct families of 3-D MHD simulations of accreting gas 
near rotating BHs that produce relativistic outflows. The limited flux simulations (LFS 
hereafter) evolve from the accretion of weak dipolar loo ps of magnetic flux of the same 



orientation from a fin i te torus of gas in the init i al state f|McKinnev and Blandford 1 12009 



Beckwith et all l2008bl ; lHawley and Krolik I l2006t iKrolik et all 120051 ). In these simulations, 
only the leading edge of the poloidal field accretes, so all the accreting large scale magnetic 
flux is of one sign. In the LFS, it is the accretion rate not the initial poloidal field strength 
that regulates the long term, large scale magnetic field strength near the black hole (and 
therefore the jet power) through ram pressure in the inner disk (private c ommunication 
McKinney (2011) and Tchekhovskoy in lMartnez-Sansigre and Rawlingsl ( 120111 ))). Physically, 
this situation might be considered a brief event of like sign large scale flux that accretes to 
the BH and is maintained near the BH by the dynamics of the accretion flow that is driven 
by a persistent MRI (magneto-rotational instability). The physical appeal of LFS is that 
memory of the initial state is erased and one finds a strong correlation between jet power 
and accretion rate in accord with the findings of PR13 that were noted above. 

The second type of simulation is based on initial conditions that creat e MCAFs (mag- 
netically choked accretion) and M ADs (magnetically arrested accretion) (JMcKinney et al 
20121 ; iTchekhovskoy et alll201ll . 120121 ). A new topology emerges, islands of large scale mag- 
netic flux perforate the disk, arresting the flow and suppressing the MRI induced dissipation 
in these regions. Thus, one would expect the luminosity to be suppressed compared to 
the LFS q It was shown in PR 13 that ~ 24 hours before an ejection Lboi is low and 



1 Simulations of MADs that were reported in lPunslv et all (|2009i ) are generally subsonic and do not produce 
significant gas heating from shocks. 
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jetted emission is minimal, thus it is not an MCAF/MAD state (which have maximal jet 
efficiency). If an MCAF/MAD switches on to drive a major flare, it would imply that 
the switch-on of an MCAF/MAD jet is preferentially associated with an increase of X-ray 
luminosity to near the highest historic (non-transient) levels. But this circumstance is con- 
tradicted by the implications of our numerical simulations th at radiative efficiency s h ould b e 
suppressed in MCAF/MADs (suppression was also noted in ISikora and Begelmanl ( 120131 ) ). 
Since MCAF/MADs simulations do not appear to be representative of the high luminosity 
during major ejections, in this article we concentrate on the LFS that can be consistent with 
the dynamics of GRS1915+105. 

In these simulations, the jet power, Q, is expressible in terms of accretion rate onto 
the BH, Mc 2 . Thus, one can compare different simulations and one can compare to the 
observations if the radiative efficiency of the accretion flow due to viscous dissipation r)th, 
(Lboi = VthMc 2 ) is known. In the numerical models, radiation effects are simulated by ad 
hoc cooling functions that are based on the local turbulent dissipation driven by MRI. In 
spite of a physically incomplete methodolog y for treating the e missivity of the accretin g 
gas, the 3-D simu lations of accretion disks in iPenna et all ( 120101 ) ; iNoble and Krolikl ( 120091 ) ; 
Noble et all (120 111 ) have been used to estimate a value, 



1. 07/ N t < ??sim = rjth < L277NT 



(5) 



wher e the range of results are expressed in units of 7/nt 5 the value from (JNovikov and Thorne 
19731 ). The disk thickness in these simulations is in the range 0.05 < H/R < 0.2. We 
consider this range of ?7 S i m from equation (5) as a reasonable estimate for the rjth of the 
LFS in Figure 5 since they have a compatible range of dis k thickness, 0.15 < H/R < 0.2 



(JMcKinney and Blandford 112009; lHawlev and Krolik 120061 ) . Furthermore, it was demon 



strated in the simulations of ISchnittman et all ( 20121 ) that even if the detailed radiative 



transfer results in the preponderance of disk luminosity being created by Compton scatter- 
ing in a disk corona instead of thermal emission from the disk proper, r] sim is unchanged 
to first order and is still consistent with Equation (5). This is an important detail because 
X-ray spectra during the ejection of major flares do not exist, so the relative contributions of 
coronal and blackbody components are unknown. The simulations neglect radiation p r essure 
in the disk, so rj t h in Equation (5) might not be accurate (but see ISzuszkiewicz et al I (119961 ) 
who calculate < 10% change in r)th for the relevant Eddington rates). 

Figure 5 compares the simulated data to the relationship depicted in Figure 4. The 
simulated data in the plot was presented and described elsewhere (see iPunslyl (120111 ) for 
deta ils). The only change is that the data is normalized by Lboi = VsimMc 2 instead of Mc 2 . 



m 



T he iMcKinnev and Blandf ord 



(120091 ) a/M=0.92 simulation shows an event horizon jet as 
Blandford and Znajeki (j 19771 ). They note that Q w O.OlMc 2 is similar to the 2-D solu- 



- 11 - 



tions reported in iMcKinneyl ( 120051 ). Thus, the spin dependent Q is given by equation (3) of 
McKinneyl (120051 ) . This normalized jet power is plotted in green in Figure 5. The other event 
horizon jet data (col l ective ly referred to a s Even t Horizon Jet (HK) in Figure 5) comes from 
Hawlev and Krolik I ( 20061 ); iKrolik et all ( 20051) except for the raw data of iBeckwith et al 
J2008a|) for a/M = 0.99 and a/M=0.998 JPunslvi IpOllh. The ergosph e ric di s k jet data is 
from the simulations, KDE, KDH and KDJ from lHawley and Krolik I (120061); iKrolik et al 
(120051 ) . The nature and strength of these jets was described in detail in IPunsly et all (120091 ) 
and references therein. In order to remove the artificial variation induced by computational 
grids covering differing amounts of the ergospheric volume in different simulations, a theo- 
retical fit to the data, "e r gosph eric disk with normalized inner boundary," (the red curve) 
was calculated in iPunslyl ( 1201 11 ). The fit was shown to exceed the KDE value because the 
inner boundary of the computational grid is farther from the event horizon in relative units 
(event horizon radius) than the other simulations. Thus it does not sample the entire ergo- 
sphere (which is critical for high BH spin in Boyer-Lindquist coordinates). The blue band 
in Figure 5 represents < Q > / < Lboi > from Figure 4 corresponding to the viable range 
of D that is compatible with the scaling to FR II quasars in Figure 3, 10.7 kpc < D < 11.0 
kpc (the dominant source of uncertainty is D). Agreement of the simulated data and obser- 
vations is achieved for an ergospheric disk driven jet with a black hole spin, a/M = 0.99 or 
a/M = 0.998. The theoretical fit in red continuously samples a/M below 0.99 and indicates 
agreement for a/M > 0.984. This result agrees with the value of a /M = 0.99 ± .01 that 



was estimated fr om the study of X-ray spectra of GRS 1915+105 ( JMcClintock et all 12006 
Blum et all 12009] ). 



Conclusion 



In this article, the long term and episodic behaviors of Lboi and the power of superluminal 
ejections in GRS 1915+105 were compared and contrasted with FR II quasars under the 
assumption of scale invariance of BH accretion systems. The results of Section 3 indicate 
that re-scaling L ho \ of GRS 1915+105 to a typical range of L bol of FR II quasars yields a 
consistent distribution of Q only if D > 10.7 kpc - otherwise the total energy emitted in 
the superluminal major ejections is too small. Physically, this constraint on the distance is 
equivalent kinematically to a Doppler factor < 0.48 and a bulk Lorentz factor > 3.5 for the 
approaching plasmoids. Comparison of the GRS 1915+105 data with the results of different 
3-D numerical simulations of the BH accretion indicate that GRS 1915+105 is compatible 
with both numerical models and re-scaling to FR II quasars if the numerical models contain 
an ergospheric disk jet and if a/M > 0.984, in agreement with observational results obtained 
by different groups and methods. If the analogy with FR II quasars is robust and scale 
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invariance is relevant to astrophysical BHs, the results obtained for GRS 1915+105 may 
imply that typical radio loud quasars also harbor rapidly spinning BHs. 
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The Ratio of the Power Required to Eject Plasmoid to the Accretion 




Flow Luminosity (<Q>/<L bo! >) 
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Fig. 4. — < Q > J < L ho i > as a function of the assumed distance, D, to GRS 1915+105. 
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Fig. 5. — The logarithm of "jet power/ Lboi" in the numerical simulations compared to 
log[< Q > J < Lboi >] for major flares in GRS 1915+105 based on Figure 4. The limited 
range of 10.7 kpc < D < 11 kpc indicated by the blue band, is the range of distance for 
which GRS 1915+105 can be considered a scaled down FR II quasar. The error bars and the 
range of uncertainty on the curves (the dashed curves) are based on the spread in Equation 
(5) 



